Inter-surface interaction via phonon in three-dimensional topological insulator 
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In topological insulators, the phonon plays the important role in electron-electron interaction. 
We consider the surface electric state on the three-dimensional topological insulator with inter- 
surface interaction via phonon theoretically. In this work, we derive the effective two-dimensional 
model in which two surface states of the topological insulator are coupled by the Fermi interaction. 
In the resultant model, the sign of the coupling constant changes depending on the temperature 
and the ratio of the speeds for the phonon and the surface state of the topological insulator. As 
a consequence, we show the prescription to tune the sign and the magnitude of the inter-surfacc 
interaction via phonon in a three-dimensional topological insulator. 



I. INTRODUCTION 



The three-dimensional(3D) topological insulator(TI) 
has two-dimensional(2D) massless electric states on its 
surface^— . The massless surface state is represented by 
a Dirac Hamiltonian in high-energy physics and has a lin- 
ear energy dispersion so called as a Dirac cone topologi- 
cally protected. Theoretically, without the interactions, a 
mass term dose not appear unless the system looses time 
reversal symmetry(TRS)^— . On the other hand, with 
the interactions, there is a well known dynamical mass 
generation for the massless fermion represented by the 
Dirac Hamiltonian in association with the spontaneously 
symmetry breaking (SSB) of chiral U(l) symmetr y 12 ; 13 . 
The dynamical mass appears when two Dirac cones are 
coupled attractively. In a 3D TI, the inter-surface in- 
teraction plays a important role in the dynamical mass 
generation because there is a single Dirac cone at one 
surface of the TI. Thus it is important to study the sign 
of the coupling constant for the inter-surface interaction. 

The TIs with the inter-surface interaction via elec- 
tromagnetic field have been studied theoretically in the 
sample of a thin filmi^— . Since the Coulomb interac- 
tion is always repulsive, the large gate voltage, i.e. the 
difference of Fermi energies for two surfaces, is neces- 
sary for SSB. However, the interaction via phonon is 
more responsible in TIs than the Coulomb interaction be- 
cause the Coulomb interaction is suppressed in TIs with 
a much larger dielectric constant than that in the vac- 
uum. The effects of the electron-electron interaction via 
phonon have been one of the hot issues in the study of 
3D TIs2ir— . However, the effects of the inter-surfaces 
interaction via phonon between two surfaces are lacked 
in the study of phonon in a TI. 

In this paper, we discuss the TI with the inter-surface 
interaction via phonon and the condition for the inter- 
surface interaction to be attractive. We construct a 
simple 2D model from the two surfaces interacting via 
phonon in a 3D TI and calculate the coupling constant 
with perturbation expansion. As a result, we find the 
inter-surface interaction via the phonon change to be at- 
tractive if the phonon speed comes close to the Fermi 
speed of the surface state. 



This paper is organized as follows. In Sec. II, we ex- 
plain our theoretical model. The inter-surface interac- 
tion is calculated numerically in Sec. III. The conclusion 
is given in Sec. IV. 



II. THEORETICAL MODEL 



We consider two pseudo (2+l)-dimensional surface 
states which are staying at the opposite surface of a 3D 
TI shown in Fig. Q] The surface state is constructed of a 
two-dimensional wave component yb = (£, if) and a wave 
component damping in the bulk of the TI. We assume 




FIG. 1. The schematic of two surface states is shown. The 
£ and r\ which are two-dimensional component of the surface 
states in the plane sample of TI. Two cones represent the en- 
ergy dispersions of the surface states. The arrows on the cones 
show the direction of spin direction locked by the momentum. 



that the penetration length A is a constant value for each 
surface states and that the surface state ^ is represented 



in thrcc-dimcnsion by, 
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where £ and 77 are the two-dimensional wave components 
of two-surface states shown in Fig. [1] In the vicinity 
of the Dirac point in two-dimensional momentum space 
parallel to the surface, the Lagrangian of the surface state 
is represented as Wcyl field by using a notation of ^> = 



7 



<*Ce 



i*7^9u* 



(4) 



where the index which appears twice in a single term 
means the summation for /j, = 0, 1, 2, 7^ is Dirac gamma 
matrices for // = — 3, and 5, <9 M represents the derivative 
for time <9o = <9o = d/dt and in-plane axis dj = vfOj = 
vpd/dxj for j = 1 and 2 with the Fermi velocity up. 
Two surface states £ and 77 are separated in a^-axis of 
a distance L. We consider a thick enough TI so that 
X/L <C 1 is satisfied and ignore a direct hopping between 
the surface states. Thus the coupling between the surface 
states is dominated by a phonon field. 



A. Effective 2D Lagrangian of electron-phonon 
coupling 

Let us derive a 2D effective electron-phonon interaction 
Lagrangian from three-dimensional model where the two 
surface states are interacting by the three-dimensional 
phonon. In the bulk of the TI, we assume the electron- 
phonon interaction can be ignored when the Fermi energy 
lies in bulk gap of the TI. The Lagrangian density of the 
phonon can be represented by use of the scalar boson 
field (f> and it is the Klein-Gordon Lagrangian density as, 

if p = (9o$) 2 - £ (v B d^) 2 -^ 2 , (5) 
3=1-3 

where vb and m are the speed and the mass of boson 
respectively. The mass appears for an optical phonon or 
an acoustical phonon under the pressure^. The free bo- 
son has the eigen energy uik = \/vB 2 \k\ 2 + m 2 with the 
momentum k. In the point of view for a lattice vibration, 
the displacement vector u is obtained by the boson field 
e.g., u = V0 for an acoustical phonon. The scalar field 
$ can be represented by a product of the wave compo- 
nent along X3-axis and the in-plane wave component of 
<j>(t,xi,X2,ka) as 



*(*,*) 



dk 3 



cj>(t,xi,X2,k 3 )—=smk 3 xs, (6) 



under the Dirichlct boundary condition of $(£, x) = at 
£3 = 0. We discuss the coupling between the Weyl field 
and the phonon field in the coordinate system with the 
boundary at x 3 = 0. 

The Lagrangian of the interaction contains the overlap 
between the wave functions of the Weyl field and the 
phonon field along a^-axis. 



(7) 



L ep = I d 3 xU{^(x)ax)f+ 2 (x 3 ) 

+i 1 \x)i 1 {x)fJ(x 3 )}<S> 1 



where the integral is defined in (#1,2:2, £3) and (,(x) = 
C(t, xi, X2) for ( = £ and 77. In the following, we omit the 
valuable x in the integrand when the valuable coincides 
with the integral valuable for convenience. The overlap 
of the wave functions of the Weyl field and the phonon 
field along a^-axis can be calculated as 



X = „ — C / dx 3 sin ^3X3 exp 
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under the limit of a small penetration length A as L/X — > 
00 and 1 + k 3 X — > 1. The interaction Lagrangian L ep is 
represented by 



L ep = / dxidxz 4>{g£, + 77+77}, (9) 

4>(t,xi,x 2 ) = dk 3 x(h)U(k 3 )(j)(t,x 1 ,X2,k 3 ). (10) 



We obtain the effective Lagrangian in 2D space (xi,X2)- 
In the following, the integral is defined in x = (t,x\,X2) 
as d 3 x = dtdx\dx2- 



B. Inter-surface interaction of Weyl field 

In this paper, we focus on representing the inter- 
surface interaction as 



V 



-9 £ f £(z) rfr)(x), 



(11) 



without the explicit phonon momentum and energy. As 
a starting point, we consider the Lagrangian of the inter- 
surface interaction intermediated by a phonon repre- 
sented by 



~&'u 



dk 3 u 2 ^i{x)D(x - x')CH(x'), (12) 



with u = xU and a phonon propagator D(x — x') (see 
Appendix [A"|) . The Lagrangian of interaction can be in- 
terpreted easily as the exchanging the phonon and shown 
as a Fcynman diagram in Fig.[5]-(b). 
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FIG. 2. The Feynman diagrams of (a) the vertex function 
in Eq.© and (b) the lowest order term for electron-electron 
interaction via the phonon exchanging are shown. The label 
£ denotes the surface Weyl field for f = £, rj. In Fig. (b), £ 
and f represent opposite surface field operator, i.e. (£, £) = 
(£, 77) or (£, C) = (?7, £)■ This diagram represent the first-order 
effective inter-surface interaction through a phonon. 



Generally, the effective action T is summation of 
the Lagrangian in Eq. ((U), (JT3J) and one-particle irre- 
ducible(lPI) vertex functions kI 1pi 1 which contains the 
loop diagrams 



1* X I <=>£• g J oi j j" 
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(13) 



In this paper, we employ the effective potential instead 
of the one-loop 1PI vertex function shown in Fig. [31 The 
effective potential is 



V 



-^lE 



,[ipi] 



(14) 



under the condition in which the initial sate and the final 
state have a same energy and momentum i.e., we assume 
the coupling constant of interaction between any states 
£ and rj in momentum space to be a single coupling con- 
stant g as 

V = -s£(pi)f(pi +p')v(P2)v(P2 -p') (15) 

where p = (po , p) are the energy po and the momentum 
p. The coupling constant g is calculated with p' = 0. 
The assumption is applicable in low energy theory when 
the interaction contains no divergent term in the infrared 
region. The effective potential in Eq. (|15[) is transformed 
to the inter-surface coupling in Eq. dill) by Fourier trans- 
formation. 

To calculate the effective potential, we introduce the 
propagators of the Weyl field and those of the phonon 
in order to calculate the effective potential in momen- 
tum space. The propagator of the Weyl field G(po,p) 
and the phonon D(ko, fe) in energy-momentum space are 
represented by 
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FIG. 3. The Feynman diagrams of the second order inter- 
surface interaction are shown. 



where a v is 2 x 2 unit matrix for v — 0, Pauli matrices 
for v = 1 — 3, pa and p = (pi,P2) are energy and two- 
dimensional momentum. 

In this work, we consider two types of coupling between 
a phonon and an electron, the deformation and polar 
coupling 2 -. Another type of coupling called as the piezo- 
electric coupling is suppressed by inversion symmetry of 
the material s 22 ' 29 and, ignored in this paper. The two 
couplings are represented by the different vertex func- 
tions Mp and Mp which are the Fourier component 
of u as 



(Mf)) 2 =C a 

(Mf)) 2 =c 2X V«B 2 N 2 + 



2 vp\p\_ 

PD 



(19) 
(20) 



where C\ and C2 are the coupling constants and the su- 
perscripts (D) and (P) denote the deformation and polar 
couplings respectively. The coupling constants are rep- 
resented by 



Ox =^R 2 

v F 

I — -- 



(21) 
(22) 



where R is the deformation constant, pn is a Debye mo- 
mentum and £oo, and e are the dielectric constants for a 
high-frequency and a low-frequency, respectively^. We 
consider both couplings in this work because it is unclear 
which coupling is more dominant in TIs. The deforma- 
tion coupling is responsible for the interaction between 
the electron and the acoustical phonon. On the other 
hand, the polar coupling accounts for the interaction be- 
tween the electron and the optical phonon. The some 
experiments suggest that the optical phonon play a dom- 
inant rol o 30 ' 31 , however the other experiment^ 2 , and the 
theory of the acoustical phonon^have a good agreement 
in the life-time of a quasi-particlc in the surface of TIs. 



C. Deformation coupling 

At first, we consider the deformation coupling. In no- 
loop calculation, the deformation coupling is vanished as 

^mt =S (0) ^b)C t C(-p) (23) 

.9 (0) = - lim / dk 3 M^ 'D( Po ,p) = 0, (24) 

where the superscript of g^ denotes the number of loops. 
Thus we must consider the one-loop effective potential 
V} and V^ in which the subscript denotes two labels in 
Fig. OH The one-loop effective potentials are represented 
by 
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(26) 



where £ and r\ are independent of (po,p). The whole 
effective potential can be represented by 



V< a 



dk 3 



Mpl'itfTtVtftTpT, 



(V^) 3 P 

+ !H 2 r2(p)?a i tT t i<r i Ti), (27) 



where the terms with a odd number of momentums pi are 
ignored here because they will be vanished at the momen- 
tum integration. Using the Matsubara method, we can 
calculate the effective potential in a finite temperature 
T = P' 1 as 



V M = -gM?t rjtf,, 



(28) 



where the second term in Eq. (|27p is cancelled out. The 
effective coupling constant g^ 1 ' is 

<?« =\ J 'dk 3 I 'dppMM^ffTiip) (29) 

ri(p) =r FB ( P ) + r F ( P ) + r B ( P ) 
r FB (p)=p[D A (p) 2 + D R (p) 2 } 

x [n F (e p )n F (-e p ) +jib(w p )jib(-w p )] (30) 
T F {p) =- — ([1- 4e p 2 D A {p)}D A (p) 2 n F (e p ) 



-[I - 4e p 2 D R (p)]D R (p) 2 n F (-e p )j (31) 

r B (p) =— ([1 + 4 Wp 2 ^ j4 (p)]^ j4 (p) 2 n B K) 



-[I + 4VZ? fl (p)]^ i? (p) 2 7i s (- Wp ) (32) 



where w p = ^/^-B 2 !^! 2 + to2 ari d e p = v_f|p| are the ener- 
gies of phonon and surface Weyl field respectively. The 



integral range of momentum is < vp 2 + k 3 < po with 
a Debye momentum p£>. The Bose and Fermi distribu- 
tion functions tib(z) and n F {z) are represented by 



n B (z) 
n F (z) 



1 



et 3z -1 
1 

"e^ + l' 



The propagator D A / R \ (p) is defined as 



D 



A(R) 



(p) 
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(33) 

(34) 

(35) 



We calculate g*- 1 -* with an unit of £dcf = Ci Pd/v f nu- 
merically by changing three parameters /3, 7 = v B /v F 
and to under the normalization by the Debye momen- 
tum p£> and the energy ui B = *\/to 2 + v b 2 pd 2 in Fig. 2] 
In the materials of TIs, 7 takes a small value of 10 -2 ap- 
proximately i.e., the phonon speed is much smaller than 
the surface fermion speed v B <C Vp. For small 7, the 
phonon contributes to the repulsive inter-surface inter- 
action shown in Fig. HJ-(a). Increasing the speed ratio 
7, however, the region of the repulsive inter-surface in- 
teraction is decreasing in Fig. HJ-(b). Eventually, the at- 
tractive inter-surface interaction becomes dominant un- 
der any mass to and temperature /3 _1 in Fig. HJ-(c). 



D. Polar coupling 

Next, we consider the polar coupling and obtain the 
finite attractive interaction in the no-loop calculation as 



ff <°> = lim fdk s M^ 2 D(p ,p) = 2C 2 «, 



(36) 



8v B s 



(n 2 2^ 4, V B p D +U} D 
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where n is always positive value. The no-loop term of 
the phonon has opposite sign to that of the photon(sce 
Appendix [Bj . We also calculate the one-loop correction 
gt 1 ' with an unit of e po i = C2 Pd/vf in this coupling 
shown in Fig. [5j For the polar coupling, the crossover 
between a repulsive coupling constant and a attractive 
coupling constant occurs at a lower 7 than the deforma- 
tion coupling. 

We should point out that it is difficult to decide which 
interaction via deformation or polar coupling is more 
dominant because the coupling constants Ci in Eq. (|I9[) 
and (|20[) arc unknown in TIs. However, these results 
with two couplings show that there is the attractive inter- 
surface interaction when the system has the slow Dirac 
mode or the fast phonon mode i.e., the speed ratio 7 
is larger than the ordinary 3D TIs of Bi2Sc3, Bi2Te3 
and Sb2Te3. In Figure El the coupling constant g^ is 
shown as a function of 7 and the temperature T = /3 _1 . 
The coupling constant g 1 - 1 ' must become attractive if the 
speed ratio 7 is larger than 0.25. 
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FIG. 4. The coupling constant g"' with an unit of Sdef = 
Ci 2 pd/vf of the inter-surface interaction via deformation in 
Eq. [28] are shown. The potential becomes attractive in red 
region and repulsive in the blue region and the zero-potential 
is represented by the green line. 



FIG. 5. The coupling constant g^ 1 ' with an unit of e po i = 
C2 2 Pd/vf of the inter-surface interaction via the polar cou- 
pling in Eq. [28] are shown. The potential becomes attractive 
in red region and repulsive in the blue region and the zero- 
potential is represented by the green line. 



At the end of this section, we show the resultant La- 
grangian density describing the opposite surface is repre- 
sented by 



££ 



-i^d^x) + g{ipijj{x) + ^7 5 V>(z))- (37) 



This model is purely two-dimensional x = (t,xi,X2) and 
the interaction is a summation as g = g(°> + g 1 - 1 ' . The 
model Lagrangian density is so called as the Nambu- Jona 
Lasino model which represent a dynamical mass genera- 
tion for a large positive g in association with the spon- 
taneously symmetry breaking of chiral C/(l) symmetry 
breakin g 12 ' 13 . We have shown that g becomes positive 
with a larger 7 than ordinal TIs. Therefore, the intrinsic 
dynamical mass generation via a phonon is possible to be 
appeared in the TI when the TI has a slow Dirac mode 



on its surface or a fast phonon mode in the bulk of TI. 



III. CONCLUSION 

We consider the inter-surface interaction via the 
phonon which propagates in the bulk of 3D TI. When 
the phonon speed is much smaller than the Fermi speed 
of the surface Dirac particle, the phonon result in a repul- 
sive inter-surface interaction. On the other hand, when 
the phonon speed catch up to the Fermi speed, the all 
phonon modes enhance an attractive inter-surface inter- 
action. The resultant effective Lagrangian can be repre- 
sented by the Nambu- Jona Lasino model with a attrac- 
tive interaction. 
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FIG. 6. The coupling constant of the inter-surface interaction 
via the polar coupling in Eq. [28] are shown. These figures 
calculated by use of (a) deformation and (b) polar coupling 
with m = 10 _3 [o;_d]. The potential becomes attractive in red 
region and repulsive in the blue region, and the zero-potential 
is represented by the green line. 
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Appendix A: Effective action of Weyl field 

In this section, we derive the effective inter-surface in- 
teraction like Eq. (fT2"|) from the Lagrangian of Eq. (J5J) and 
([9]) by the path-integral method. The effective action T 
which contains the effective interaction is acquired from 
the generating functional of a propagator Z represented 
by, 

Z[z,z,J] = J @ip@i>@(t>ex-pi J d 3 x£? (Al) 

Jz? = Jz? e + Jz? p + u a (^ ip)4> + J<P + zip + zrj>, 

(A2) 

where z(x), z(x) and J(x) are the virtual external sources 
and we introduce the coupling constant uq. The ef- 
fective action is the generating functional of the vertex 
function^. To omit an explicit phonon field, we consider 
that there is no phonon source as J = and calculate 



Z[z,z] = / ^V^exp 

J2? e ' =J^ e + Z1p + Zip, 



d 3 x^f! 



Z \j(x)], (A3) 



z o[j{x)} =exp 



—i I d x I d x'j(x)D(x — x')j(x') 



j(x) =u(ipj tp)(x), 

with the phonon propagator D(x). The calculation is 
same as the case of a non-interacting scalar boson with 
external source of j(x) and provides the generating func- 
tional of phonon propagator as Zq\j(x)]. We obtain the 
four point term in Zq which represents an effective inter- 
action in Weyl field, 



*U = -J(X)D(X - X')j(x'). 



(A4) 



The interaction can be interpreted easily as the exchang- 
ing the phonon. The Lagrangian contains an intra- 
surface interaction proportional to the fourth-order term 
of £ or r] in Eq. (fTTj) and an inter-surface interaction which 
consist of the product of the quadratic term of £ and the 
quadratic term of ?/. The effective action can be derived 
as a summation of one-particle irreducible vertex func- 
tions which consist of four-field potential in Eq. (|A4[) . 



Appendix B: The Coulomb interaction 

The interaction Lagrangian via a photon in Feynman 
gauge is represented by 






-9^ 



Po< 



\P\ 



IE 



(Bl) 
(B2) 



where the speed of light is unity c = 1 and g^ v is the 
metric tensor as, 



g^ =diag[l, -1,-1,-1], 



(B3) 



with [i, v — — 3. Using the longitudinal coiuponent 
given by -Doo (p) > the interaction without energy transfer 
Pa = can be represented by 



~£\x 



(te + Jri-^tft + tfr)) 



\P\ 



w [ EcWc+aeW,, 



(B4) 

(B5) 



It is obviously for the first and second terms are the 
intra-surface and inter-surface Coulomb interactions re- 
spectively. This interaction is repulsive in contrast to 
the phonon in Eq. Q36[) because the photon propagator 
in Eq. (|B2[) is opposite sign to the phonon propagator in 

Eq. CGD. 



1 L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

2 L. Fu and C. L. Kane, Phys. Rev. B. 76, 045302 (2007). 

3 J. E. Moore and L. Balents, Phys. Rev. B. 75, 121306 
(2007). 

4 Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S. K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 
Zhang, I. R. Fisher, Z. Hussain, and Z. X. Shen, Science 
325, 178 (2009). 

5 Y. Tanaka, T. Yokoyama, and N. Nagaosa, Phys. Rev. 
Lett. 103, 107002 (2009). 

6 D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, 
J. Osterwalder, L. P. J. G. Checkelsky, N. P. Ong, A. V. 
Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. 
Cava, and M. Z. Hasan, Nature 460, 1101 (2009). 

7 Q. Liu, C.-X. Liu, C. Xu, X.-L. Qi, , and S.-C. Zhang, 
Phys. Rev. Lett. 102, 156603 (2009). 

8 Y. S. Hor, P. Roushan, H. Beidenkopf, J. Seo, D. Qu, 
J. G. Checkelsky, L. A. Wray, D. Hsieh, Y. Xia, S.-Y. Xu, 
D. Qian, M. Z. Hasan, N. P. Ong, A. Yazdani, and R. J. 
Cava, Phys. Rev. B. 81, 195203 (2012). 

9 Y. L. Chen, J.-H. Chu, J. G. Analytis, Z. K. Liu, 
K. Igarashi, H.-H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 
D. H. Lu, M. Hashimoto, T. Sasagawa, S. C. Zhang, I. R. 
Fisher, Z. Hussain, and Z. X. Shen, Science 329, 659 
(2010). 

10 L. A. Wray, S.-Y. Xu, Y. Xia, D. Hsieh, A. V. Fedorov, 
Y. S. Hor, R. J. Cava, A. Bansil, H. Lin, and M. Z. Hasan, 
Nat. Phys. 7, 32 (2010). 

11 T. Habe and Y. Asano, Phys. Rev. B. 85, 195325 (2012). 

12 Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 
(1961). 

13 Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 
(1961). 

14 B. Seradjeh, J. E. Moore, and M. Franz, Phys. Rev. Lett. 
103, 066402 (2009). 

15 J. Wang, A. M. DaSilva, C.-Z. Chang, K. He, J. K. Jain, 
N. Samarth, X.-C. Ma, Q.-K. Xue, and M. H. W. Chan, 
Phys. Rev. B 83, 245438 (2011). 



16 M. Liu, C.-Z. Chang, Z. Zhang, Y. Zhang, W. Ruan, K. He, 
L.-l. Wang, X. Chen, J.-F. Jia, S.-C. Zhang, Q.-K. Xue, 
X. Ma, and Y. Wang, Phys. Rev. B 83, 165440 (2011). 

17 G. Y. Cho and J. E. Moore, Phys. Rev. B 84, 165101 
(2011). 

18 E. G. Moon and C. Xu, Eur. Phys. Lett. 97, 66008 (2012). 

19 I. Sodemann, D. A. Pesin, and A. H. MacDonald, Phys. 
Rev. B 85, 195136 (2012). 

20 D. K. Efimkin, Y. E. Lozovik, and A. A. Sokolik, Phys. 
Rev. B 86, 115436 (2012). 

21 W. Cheng and S.-F. Ren, Phys. Rev. B 83, 094301 (2011). 

22 S. Giraud and R. Egger, Phys. Rev. B 83, 245322 (2011). 

23 Z.-H. Pan, A. V. Fedorov, D. Gardner, Y. S. Lee, S. Chu, 
and T. Valla, Phys. Rev. Lett. 108, 187001 (2012). 

24 X. Zhu, L. Santos, R. Sankar, S. Chikara, C. Howard, F. C. 
Chou, C. Chamon, and M. El-Batanouny, Phys. Rev. 
Lett. 107, 186102 (2011). 

25 S. Giraud, A. Kundu, and R. Egger, Phys. Rev. B 85, 
035441 (2012). 

26 X. Zhu, L. Santos, C. Howard, R. Sankar, F. C. Chou, 
C. Chamon, and M. El-Batanouny, Phys. Rev. Lett. 
108, 185501 (2012). 

27 O. Viyuela, A. Rivas, and M. A. Martin-Delgado, Phys. 
Rev. B 86, 155140 (2012). 

28 G. D. Mahan, Many-Particle Physics THIRD EDITION 
(Plenum, New York, 2000). 

29 B.-L. Huang and M. Kaviany, Phys. Rev. B 77, 125209 
(2008). 

K. Shahil, M. Hossain, D. Teweldebrhan, and A. Balandin, 
Appl. Phys. Lett. 96, 153105 (2010). 

31 J. Qi, X. Chen, W. Yu, P. Cadden-Zimansky, D. Smirnov, 
N. H. Tolk, I. Miotkowski, H. Cao, Y. Chen, Y. Wu, 
S. Qiao, and Z. Jiang, Appl. Phys. Lett. 97, 182102 (2010). 

32 R. C. Hatch, M. Bianchi, D. Guan, S. Bao, J. Mi, B. B. 
Iversen, L. Nilsson, L. Hornekaer, and P. Hofmann, Phys. 
Rev. B. 83, 241303 (2011). 

33 J. W. Negele and H. Orland, QUANTUM MANY- 
PARTICLE SYSTEMS (Westview Press, 1998). 



